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We examine cosmological, astrophysical and collider constraints on thermal dark matter 
(DM) with mass mx in the range ~ 1 MeV — 10 GeV. Cosmic microwave background 
(CMB) observations, which severely constrain light symmetric DM, can be evaded if the 
DM relic density is sufficiently asymmetric. CMB constraints require the present anti-DM 
to DM ratio to be less than - 2 x 1CT 6 (1CT 1 ) for DM mass m x = 1 MeV (10 GeV) with 
ionizing efficiency factor / ~ 1. We determine the minimum annihilation cross section for 
achieving these asymmetries subject to the relic density constraint; these cross sections are 
larger than the usual thermal annihilation cross section. On account of collider constraints, 
such annihilation cross sections can only be obtained by invoking light mediators. These light 
mediators can give rise to significant DM self-interactions, and we derive a lower bound on the 
mediator mass from elliptical DM halo shape constraints. We find that halo shapes require 
a mediator with mass m^>4x 10 -2 MeV (40 MeV) for m x =1 MeV (10 GeV). We map 
all of these constraints to the parameter space of DM-electron and DM-nucleon scattering 
cross sections for direct detection. For DM-electron scattering, a significant fraction of the 
parameter space is already ruled out by beam-dump and supernova cooling constraints. 
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I. INTRODUCTION 



Studies of dark matter (DM) have historically focused on particles with weak scale mass ~ 100 GeV 
The reason is not only the focus of the high energy physics community on weak scale phenomena, but also 
because the annihilation cross section for a Weakly Interacting Massive Particle (WIMP) naturally gives rise 
to the observed cold DM relic abundance. This is the so-called "WIMP miracle." 

More recently there has been a broader interest in light DM, with mass mx $10 GeV. Part of the 
reason for this interest is phenomenological. Direct detection results from DAMA [4[, CoGeNT [5, 6], and 
CRESST 7] claim event excesses that can be interpreted as nuclear scattering of DM with mass ~ 10 GeV 
(although the mutual consistency of these results is disputed). Meanwhile dark matter with masses of MeV 
has been studied as a possible explanation of the INTEGRAL 511 keV line @-14|. 

There is also a theoretical motivation for light DM, as DM with mass mx < 10 GeV appears in certain 
classes of models naturally. In supersymmetric hidden sector models, for exam ple, gauge interactions generate 
light DM masses and give rise to the correct annihilation cross section 12} [la, |l6|. The asymmetric DM 
(ADM) scenario, where the DM particle X carries a chemical potential, analogous to the baryons, provides 
another approach to light DM (see e.g. 17-2(| and references therein). In these scenarios, both DM (X) and 
anti-DM (X) particles may populate the thermal bath in the early Universe; however, the present number 
density is determined not only by the annihilation cross section, but also by the DM number asymmetry 
•qx ■ Depending on the value for r\x , the DM mass can be as low as ~ keV in ADM models [5l[ , though the 
natural scale for ADM is set by (f2cDM/^b)"T-p ~ 5 GeV. 

The purpose of this paper is to explore model-independent constraints and predictions for the asymmetric 
and symmetric limits of light DM with mass ~ 1 MeV— 10 GeV. 1 Although both phenomenological and 
theoretical considerations have motivated the study of light DM candidates, there are still a number of 
important constraints that should be taken into account in realistic model building. In general, light thermal 
DM faces two challenges: one is to evade bounds on energy injection around redshifts z ~ 100 — 1000 
coming from observations of the CMB; the other is to achieve the required annihilation cross section without 
conflicting with collider physics constraints. 

CMB data from WMAP7 strongly limits DM annihilation during the epoch of recombination, and excludes 
symmetric thermal light DM with mass below ~ 1 — 10 GeV if the annihilation is through s-wave processes 



22H24j . The CMB bounds may be evaded in the symmetric case if DM dominantly annihilates to neutrinos or 



if its annihilation is p-wave suppressed. When the DM relic density is asymmetric, DM annihilation during 
recombination can be highly suppressed if the symmetric component is sufficiently depleted, providing a 
natural way to resolve the tension from CMB constraints for light DM scenarios. Unlike the case of symmetric 
DM, the CMB places a lower bound on the annihilation cross section for ADM from the requirement of 
sufficient depletion of the symmetric component. We calculate the minimum annihilation cross section 
required in order to evade the CMB bound and achieve the correct relic density simultaneously. 

However, it is difficult to achieve the needed annihilation rate to Standard Model (SM) particles through 
a weak-scale mediator. Null results from mono-jet plus missing energy searches at the Tevatron 25j-|27| 
and the LHC [H Hjj strongly constrain such a mediator if DM couples to quarks and gluons. Meanwhile, 
the moncHphoton plus missing energy search at LEP sets limits on the coupling between DM and charged 
leptons [3Cj via such a heavy state. These collider constraints are so strong that the annihilation through 
an off-shell heavy mediator is generally insufficient for ADM to achieve the correct relic density and evade 
the CMB constraint, if the DM mass is below a few GeV. One way to evade the collider constraints is to 
invoke a light mediator with mass much less than ~ 100 GeV. In this case, DM can annihilate to SM states 
efficiently via the light state without conflicting with collider bounds. Furthermore, if the mediator is lighter 
than the DM, a new annihilation channel opens and DM can annihilate dominantly to the mediator directly. 
In this limit, the mediator particle may couple to the SM sector rather weakly. 



1 For DM much lighter than ~1 MeV, DM can only annihilate to neutrinos, new light states that remain relativistic through 
matter-radiation equality, or hidden sector forces that decay invisibly. In this case, the CMB and collider bounds discussed 
here do not apply. 
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The presence of the light mediator has various implications for DM dynamics in galaxies and for cosmology. 
The light mediator may give rise to significant DM self-interactions (i.e., DM-DM scattering); this is true 
in both the symmetric and asymmetric limits, since the light state mediates DM-DM interactions as well 
as anti-DM and DM interactions. These interactions leave footprints in the DM halo dynamics. There are 
limits on the DM self-interaction cross section coming from observations of the Bullet Cluster, elliptical 
galaxy clusters, and elliptical DM halos. We combine these with the relic density constraint to place a lower 
bound on the mediator mass ~4x 10~ 2 MeV — 40 MeV for DM masses in the range ~ 1 MeV — 10 GeV. 
We assume this massive mediator decays to SM relativistic degrees of freedom in the early universe to avoid 
the overclosure problem, and derive conditions for thermalization of the DM and SM sectors. 

These astrophysical and cosmological constraints can be applied to the parameter space of scattering rates 
in direct detection experiments. We consider DM-nucleon scattering for DM masses of 1 — 10 GcV and 
DM-electron scattering for DM masses 1 MeV — 1 GeV. In the case of electron scattering, we combine 
the astrophysical and cosmological constraints with bounds from beam dump experiments and supernova 
cooling, which exclude a large region of the allowed parameter space. In addition, the predictions are very 
different dependent on whether the mediator is heavier or lighter than the DM. 

The rest of paper is organized as follows. In Section II, we present the relic density calculation for DM 
in the presence of a chemical potential. In Section III, we study the CMB constraint on ADM models and 
derive the annihilation cross section required to evade the CMB bound. In Section IV, we examine current 
collider physics constraints on the DM annihilation cross section. In Section V, we study the elliptical halo 
shape constraint on the mediator mass. In Section VI, we map out the parameter space for DM direct 
detection. We conclude in Section VII. 



II. RELIC DENSITY FOR SYMMETRIC AND ASYMMETRIC DARK MATTER 



Our starting point is to establish that the correct relic density of £7cdm^ 2 = 0.1109 ± 0.0056 [31| can 
be obtained, where we assume that the annihilation cross section (av) and the asymmetry rjx are floating 
parameters. 

In the usual thermal WIMP scenario, the correct relic density is determined by DM annihilation 
until freeze-out. For Dirac DM in the symmetric limit, the cold DM relic density is ^cdm/i 2 ~ 
0.11 (6 x 10 _26 cm 3 /s) /(crv). DM may also carry a chemical potential which leads to an asymmetry be- 
tween the number density of DM and anti-DM. In this case, when the DM sector is thermalized, the 
present relic density is determined both by the annihilation cross section and the primordial DM asym- 
metry rjx = (nx — n x )/s, where nx, n x are the DM and anti-DM number densities and s is the entropy 
density. In the asymmetric limit, neglecting any washout or dilution effects, the correct relic density is 
obtained for a primordial asymmetry given by 

^CDM Pc , 1S 
VX ~ , (1) 

rax s 

where sq s» 2969.5 cm~ 3 and p c s» 1.0540/i 2 x 10 4 eV/cm 3 are the entropy density and critical density today. 
In the asymmetric limit, the annihilation cross section is sufficiently large that the thermally-populated 
symmetric component is a sub-dominant component of the energy density today. 

Depending on the strength of indirect constraints on DM annihilation, light DM scenarios must interpolate 
between the symmetric and asymmetric limits. We thus require precise calculations of the present anti-DM 
to DM ratio = Q, x JQ,x, which controls the size of indirect signals from DM annihilation. Note that 
is related to the absolute relic densities by 

1 Vxm x So r]xmxs 
= -j— , = 73 > K 2 -) 

± Too Pc -L ^*oo Pc 

and the total CDM relic density is SIcdm = + &x ■ 

To compute we solve the Boltzmann equations for nx , n x freezeout in the presence of a nonzero 
chemical potential [32]. In this work, we focus on the case where DM is in thermal equilibrium with the 
photon thermal bath through freezeout. In general, this assumption may not hold if there is a weakly coupled 
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light mediator coupling DM to the SM. We leave the more general case for future work 33[, noting that the 
effects on the relic density are up to 0(10), depending on mx- 

The coupled Boltzmann equations for the species n + = nx and n_ = n x are 

— — = —3Hn± — (av) [n + n_ — n^'nl 9 ] (3) 

where (av) is the thermally-averaged annihilation cross section over the X and X phase space distribu- 
tions [34]. The Hubble expansion rate is H « 1.66^/g^T 2 /M pl where M pl « 1.22 x 10 19 GeV is the Planck 
mass and g e f{ is the effective number of degrees of freedom for the energy density. If there is a primordial 
asymmetry in X number, then there is a nonzero chemical potential u which appears in the equilibrium 
distributions as n!|? = e ±fi / T n eq . Here n eq is the usual equilibrium distribution with u = 0, and thus 
n^nl 9 - (n e ?) 2 . 

We then take the standard definitions x = mx/T and Y± — n±/s, where s = (2ir 2 /45)h c s(T)T 3 is the 
entropy density and h e g(T) is the effective number of degrees of freedom for the entropy density. We write 
the annihilation cross section as (av) — aoX~ n , with n = and n = 1 for s-wave and p-wave annihilation 
processes respectively. Then simplifying Eq. (|3]) gives 

d ^t = -^V9:{Y + y—(Y^f), ( 4 ) 

where A = O.264M p vmjf0o an d Y eq — ^-^^(g/h e g)x 3 ^ 2 e^ x = ax 3 ^ 2 e~ x . The effective number of degrees of 
freedom is ^ = ^ (l + 3^^) 0. 

After being generated at some high temperature, the DM asymmetry is a conserved quantity, so we have 
the constraint 

vx = y+ - Y- (5) 

which is constant at any given epoch. 2 In order to impose this condition on our numerical solutions, we 
define the departure from equilibrium A by Y± — Y± q + A, and instead solve the (single) equation for A. 
It is helpful to present approximate analytic solutions in the limit of constant ^fgl 32, 381 3^]. Eq. (f4]) 



can be solved analytically at late times when (Y eq ) 2 becomes negligible. In this limit, using Eq. ([5]), we can 
integrate Eq. (jH) separately for X and X to obtain 

1±(00) ~ ^ ; . (6) 

1 — [1 =F Vx/Y±(x f )] e* vxX ^ x f /(n+1) 

These solutions also apply for the symmetric case in the limit of r\x —> 0. We take the freezeout temperature 
Xf = mx I Tf as derived in [38| : 

, u 1 1 \ / — \i 1 1 In 2 [(n + l)VgaA] 

x f ~ In [in + \)JgZa\ + - In „ - -. — 5-. (7) 

2 ln 2 " +4 [(n + l)^ a A] ~(^g;) 2 l(n + l)\ V x/2} 2 

Using Y±(oo) given in Eq^. ([5]), we can obtain the present ratio of the X to X number densities: 

Y-, v Y-(xf) ( -vxxyg: \ fa . 

While we can obtain a precise analytic result for r(xf) — Y—(xf)/Y+(xf), it turns out that the consequence 
of neglecting the (Y eq ) 2 in the late-time solution can almost exactly be accounted for by simply setting 
r(xf) = 1. This gives numerically accurate answers over a wide range of r\x arid (av) as discussed in [38j |. 
Note that the solution here only converges when r\x A is small enough yfglr\x^ < 2a;^ +2 . 



2 We assume there is no Majorana mass term for DM, and thus X — X oscillation [2lL l35|- [3^l does not occur. We also assume 
there is no entropy production in this case and there are no DM-number violating interactions at these temperatures. 
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III. CMB CONSTRAINTS 



For both symmetric and asymmetric thermal DM, the DM particles must have a sufficiently large annihi- 
lation cross section in order to achieve the correct relic density. This annihilation may have many indirect 
astrophysical signatures; among these, the most robust prediction (or constraint) is the effect of DM anni- 
hilation on the cosmic microwave background (CMB) |40j], since the effect only depends on the average DM 
energy density. We first summarize recent studies of CMB constraints on DM annihilation, and then discuss 
scenarios which naturally evade these constraints for light DM, focusing on the asymmetric DM scenario. 

Energy deposition from DM annihilation distorts the surface of last scattering, which affects the CMB 
anisotropies and is thus constrained by WMAP7 data. CMB constraints become increasingly severe for 
smaller DM masses: the energy released in DM annihilations scales as ~ mx(nx) 2 ~ PcDM/ m x, where 
Pcdm is the average energy density in DM. This implies the effect of DM annihilation on the CMB scales 
as ~ (av)/mx- Though the precise bound depends on the mass and annihilation channels, WMAP7 limits 
the amount of annihilation during recombination to below the thermal relic annihilation cross section if 



mx < 1 — 10 GeV [22J, |23|, |4l|, |42j . Furthermore, Planck data can improve these constraints by up to a factor 
of 10. 

For self-annihilating DM particles such as Majorana fermions or real scalars, the energy deposition rate 
per volume at redshift z is 

——(z)=p c n cDM {l + z)°f(z) , (9) 

did V mx 

where p c is the critical density at the present time, (crv) cmb is the thermally-averaged annihilation cross 
section at the epoch of recombination, and f(z) parametrizes the amount of energy absorbed by the photon- 
baryon fluid at redshift z, relative to the total energy released by DM annihilation at that redshift. 

The quantity f(z) gives the efficiency of energy deposition at redshift z and thus depends on the spectrum 
of photons, neutrinos and resulting from DM annihilation. In general, the dependence of f(z) on z is 
mild [41|, and an excellent approximation is to take f(z) = /ewiMp(z) where / is a constant and ewiMp(z) 
is a universal function for WIMP DM [24|. In addition, to leading order / ~ (1 — f v ) 23 1, where f v is the 



fraction of energy going to neutrinos per annihilation. For DM annihilation channels to charged lepton or 
pion final states, / ~ 0.2 — 1; here annihilation only to can give / ~ 1. 

There is also some mild mx dependence in f(z) (or /), since the spectrum of DM annihilation products 
depends on mx- Ref. [41[ computed detailed efficiency curves f(z) for mx > 1 — 10 GeV, depending on 
the channel. However, the observed trend is that efficiency does not depend strongly on mass in the range 
1-1000 GeV, and furthermore increases for lower mass. 3 We will extrapolate results to mx < 1 GeV; we 
expect this is a conservative approach. 

The WMAP7 limit on DM energy injection at the 95% C.L. can be written as [22] 

, Wcmb < 2.42 x IP' 27 cm 3 /s 
m x GeV 

This bound 4 as given assumes DM particles are self-annihilating, i.e. Majorana fermions or real scalars. For 
DM candidates that are Dirac fermions or complex scalars, as in ADM scenarios, the energy injection rate is 

dE . . 2 2 g ,(<Jv) C MB 

-\z) =2p c n GI}M — —— —(l + z)°f(z)— , (11) 



dtdV y ' rc OUM (l+roo) 2 m x 



3 Above mx,rn,p > 1 MeV, most of the annihilation products rapidly cascade down to lower energies and the efficiency / is 
only mildly sensitive to the initial energy spectrum of annihilation products (normalizing for the total energy). However, 
photons in the range ~ 0.1 — 1 GeV deposit their energy relatively inefficiently. For annihilation of sub-GeV scale DM, 
typically a smaller fraction of the total energy goes into photons of these energies, which increases the total efficiency slightly. 
We thank Tracy Slatyer for this point. 

4 Note: the results of |23H are slightly weaker by a factor of 1.2-2. 
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Figure 1: WMAP7 95% C.L. constraints on the DM annihilation cross section and mass for asymmetric dark matter 
and s-wave annihilation. We show constraints for various values of r = r^, = Q. x /Q.x, the anti-DM to DM ratio at 
the present time. The shaded region (blue) is excluded by the WMAP7 data, with different shades corresponding to 
different r x . Along the horizontal contours of constant r are the values of (av) where the correct relic density can 
be obtained for an efficiency factor / = f . The turnover around mx ~ 10 GeV comes from the drop in SM degrees 
of freedom when the universe has temperature ~ 1 GeV. The solid red line is the intersection of the WMAP7 and 
relic density contours: it indicates the minimum (av) needed to obtain the observed relic density and satisfy CMB 
constraints for s-wave annihilation. 



where we have used px + Px = Pcdm and = px/px- Note there is factor of 2 in the energy injection 
rate relative to the self-annihilating case, accounting for the number of possible annihilations. Comparing 
Eq. (JHJ) and Eq. (fTTj) , we can translate the bound given in Eq. (fTU)) to the Dirac fermion or complex scalar 
case: 

2roo Hqmb 2.42 x 10~ 27 cm 3 /s 
(l+r QO )2 / mx GeV 1 j 

We show this constraint for various values in Fig. ([1]); the dotted black line gives the thermal relic 
annihilation cross section in the symmetric case, where we have solved for the relic density numerically and 
taken / = 1. 

ADM can evade CMB bounds while still allowing s-wave annihilation. 5 The CMB bounds do not com- 



5 In the symmetric limit, one can evade the CMB bounds if DM annihilates via p-wave suppressed interactions. Then 
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Figure 2: (Top) Minimum (av) for efficient annihilation of the symmetric component in an ADM scenario, such that 
CMB bounds can be evaded, for two different values of the efficiency /. The black dotted line gives the thermal relic 
(av) for the symmetric case. (Bottom) The corresponding maximum allowed r^, the anti-DM to DM ratio at the 
present time. 



pletely disappear in the ADM scenario, however, because there is a small symmetric component of DM 
remaining, r^, the size of which depends on (av). Because of the exponential dependence of r m on (av), as 
shown in Eq. (j8]), the CMB constraints lead to a lower bound on (av). This is shown in Fig. (p}, where we 
map out the constraints in the (ctv)cmb an d Tnx parameter space, computing the relic density numerically 
and applying the constraint in Eq. (|10[) . The solid line (red) gives the resulting lower bound on f(av)cMB- 
This lower bound on f(av)cMB translates to an upper bound on the residual symmetric component, r^, as 
shown in Fig. @. We give analytic approximations to these numerical solutions next. 

When Too < 1, we can ignore the X contribution to the total relic density, and the DM asymmetry 
parameter rjx is set by i]x ~ ^cdmPc/ '(mxso)- For a given rjx, the required annihilation cross section at 
freezeout to achieve a particular residual symmetric component, r^, can be obtained by rewriting Eq. ([5]) 
as 

( av)f ^ £52 l n ( ± 

0.264:ilcBMPcy/g7j M pl \ r oo 

~ Cf x 5 x 1CT 26 cm 3 /s x ln (^-j , (13) 

where c/ = (|jy) ( ^ = j is an 0(\) factor. We show the numerical result as the horizontal contours of 
constant r m in Fig. ([T}; for mx < 1 GeV we obtain a good approximation to the numerical solution by 



(™)cmb — ("CMb/"/) 2 (c^)/ an( l since i>cmb ~ 10 — 8 while Vj ~ 0.3, the annihilation cross section at recombination 
is highly suppressed and WMAP constraints are substantially weakened. An increased branching ratio to neutrinos (smaller 
/) can also alleviate the tension with CMB data for light DM. 
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taking c/ = 1. On the other hand, the CMB bound on the annihilation cross section when 1 is 

, s 2.42 x 1CT 27 cm 3 /s ( m x \ f I \ 

( ™ )CMB < 27 " (icfv) [tz ) • (14) 

For s-wave annihilation, we take (av) f — (c™)cmb- Since (av) t increases with log^/r^), but the CMB 
bound on ((tw)cmb increases with I/Too, we can evade the CMB constraints by decreasing r^. For a given 
DM mass, thermal ADM is consistent with the CMB constraints if satisfies the following condition, 

\rooJ f xcf VI GcV/ 

The numerical result for this bound is shown in Fig. @; a good analytic approximation is given by < 
ro/ln(l/ro), with tq ~ 2 x 10 _2 (mjc/GeV)//. Taking / ~ 1, we can see that has to be smaller than 
5 x 1(T 3 and 2 x for m x ~ 1 GeV and 1 MeV, respectively. 

Likewise, we can combine Eq. (ITS)) and Eq. flT4"|) to place a lower bound on (<jv) /: 



h > fin (40c f fx^f) +lnln(4Q C/ / xiM) , m x < / x lOGeV. 

— 26 r-m 3 / o r ^ I r\ \ /* . . 1 nn .T7 



c f x 5 x 10 26 cm 3 /s ~ I 2 , m x > f X lOGeV. 



(16) 



Note if mjc is larger than / x 10 GeV, the CMB constraints do not apply and the annihilation cross section is 
set by the relic density requirement. The analytic approximation in Eq. (|16[) agrees well with the numerical 
results, which are shown in Fig. 

With these constraints on the minimum annihilation cross section, we now turn to discussing what classes 
of models can generate the needed annihilation cross section consistent with collider constraints. 



IV. LIGHT MEDIATORS 



Thus far, we have treated the annihilation cross section (av) as a free parameter. To proceed we must 
specify the physics that generates this cross section. First, DM may annihilate directly to SM particles 
through heavy mediators with mass greater than the weak scale. This coupling to the SM implies light 
DM can be produced in abundance in colliders. We review constraints from missing (transverse) energy 
searches at collider experiments and from direct detection experiments, which conflict with the (av) required 
to obtain the observed relic density. In this case, thermal light DM is ruled out in both the symmetric and 
asymmetric scenarios. Second, DM can annihilate via new light states which have a mass below the typical 
momentum transfer scale in the colliders. In this case, the collider constraint can be evaded. If the new 
state is lighter than DM, it can be very weakly coupled to the SM. 



A. Collider and Direct Detection Constraints on Light DM with Heavy Mediators 

In the heavy mediator case, a convenient way to parametrize the DM-SM coupling is via higher dimensional 
operators, which is valid if the mediator mass is heavier than the relevant energy scale. Here we give two 
typical examples, 



r) X ^ x f^f An XX ^ 

O x : and 2 : ^ > ( 17 ) 



where X is DM, / is a SM fermion, and A 12 are cut-off scales for 0\^- The cut-off scale, in terms of the 
parameters in the UV-complete models, is A = m^/y/gxgf, where m<j, is the mediator mass, and gx and <?/ 
are coupling constants of DM-mediator and SM-mediator interactions respectively. 

In the limit of mx 3> to/ , the DM annihilation cross sections at freezeout are given by 

N c m 2 Afc 2 -i 

H^jx^^^jz (18) 
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for C>i and O2 respectively. N'j is the color multiplicity factor of fermion /, and Xf = mx/T ps 20, with 
T the temperature. Note that the annihilation cross section through O2 is p-wave suppressed. Now we can 
estimate the limit on the cut-off scales Ai and A2 by requiring the correct relic density 



3 J VlO GeV/ V (<™) / \ x f J 

where the limit is relevant for both the asymmetric and symmetric cases. Since the annihilation cross section 
is p-wave suppressed for C 2 , we need a smaller cut-off scale to obtain the correct relic abundance. Now we 
review various constraints on the cut-off scales Ai,2- 

• Direct Detection Constraints 

If DM couples to quarks, the operators ©1,2 can lead to direct detection signals with the DM-nuclcon 
scattering cross section: a ni 2 ~ /i^/Aj 2 , and fx n is the DM-nucleon reduced mass. For a DM mass 
~ fO GeV, taking the value of Aj.,2 given in Eqs. (ITO)) and (f2"0"|). we expect the DM-nucleon scattering 
cross section to be a ni ~ 10~ 38 cm 2 and a n2 ~ 10 -36 cm 2 . However, the current upper bound on a n 



from direct detection experiments for DM with mass mx > 10 GeV is a n < 10~ 42 cm 2 43[, which is 
much smaller than the predicted values from requiring the correct thermal relic density. For DM with 
mass below a few GeV, the recoil energies are too small and direct detection bounds are currently very 
weak or nonexistent. 

Tevatron and LHC Constraints 

The DM-quark interactions given in can lead to signals of mono-jet plus missing transverse energy 
at hadron colliders, while the Tevatron data for this signal matches the SM prediction well. We require 
that Oi t 2 do not give rise to sizable contributions to this signal. The lower bounds on Ai 2 are 
~ 400 GeV and ~ 400 GeV jUSzl respectively, for DM masses mx < 10 GeV that we are interested 



in. Recent LHC results give a stronger limit on Ai > 700 GeV [29j. Therefore the Tevatron and LHC 
searches have excluded both thermal symmetric DM and ADM in the whole range of light DM if the 
DM particles annihilate to light quarks through 0\ and 02- 

• LEP Constraints 

If DM particles couples to the electron through ©1,2, the mono-photon search at LEP sets a limit on 
the cut-off scale: Ai > 480 GeV and A 2 > 440 GeV for DM mass m x < 10 GeV do|. Note the limit 
also applies to the case where DM couples to three generations of charged leptons universally. One 
may avoid the limit by coupling DM only to /i or r. However this approach usually involves model 
building complications and severe flavor constraints. 

Thus we conclude that for ©1,2, DM does not have the correct relic abundance for symmetric DM and 
ADM due to the combination of direct detection and collider constraints. The direct detection constraints 
can be relaxed by suppressing the direct detection scattering cross section; this can happen for example if the 
scattering off nuclei is velocity suppressed, notably through an axial interaction. However, the collider bounds 
are still severe for higher dimensional operators involving interactions with light quarks or electrons [25l - l30l | . 



B. Light Dark Matter with Light Mediators 



One simple way to evade the collider constraints for light DM is to invoke light mediators with masses 
much smaller than the typical transverse momentum of the colliders pt ~ O(100 GeV) (or the center of 



mass energy ~ 200 GeV for LEP) 
bounds become much weaker [27 
probed at colliders, there exists a large parameter space for light DM scenarios to achieve the correct relic 



In this limit, the effective theory approach breaks down and the collider 
29l 30l . 44 1 . In general, if the mediator mass is much less than the px 
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density. We consider a hidden sector with Dirac DM coupled to a light mediator which could be a spin-1 or 
spin-0 particle; for ease of notation we always refer to it as <p. We write the Lagrangians as 

C v = gxX^X^ +g f frf^ + m x XX + m^ M , (21) 
C s = gxXXcj) + g f ff<f> + m x XX + m 2 2 , (22) 

where is the mediator mass. We consider two cases for the mediator mass: 6 a mediator with > 2mx 
and lighter mediator with m^ < mx ■ 

In the case of px 3> vn<$> > 2mx , the DM particles can annihilate to SM particles through the s-channel 
process. There is a collider bound on gj because an on-shell mediator which decays to XX can be produced, 
potentially contributing to the mono-jet plus missing transverse energy signal. Tevatron data has been 



employed to place an upper bound on gj < 0.015/ '^/Br(<f> ~ > XX) for < 20 GeV 44|, where Br( 
XX) is the branching ratio of <f> decay to the DM pair. In this case the annihilation cross section is given by 
(<jv)v — 4a x g'jm x Nf / 'm^ and (av)s — a x g 2 fm x N'j/2m^ > Xf 1 where ax = g\j^- To see how the collider 
constraint affects the annihilation cross section in this case, we take the conservative limit gj < 0.015, setting 
Br(<j) — > XX) ~ 1. From the relic density constraint, we then obtain an upper bound on the mediator mass, 

m^lSGeV^j ( {(Tv) ) (iGevJ ■ ( 23 ) 

This bound 7 is consistent with our assumption that m$ 3> mx- 

If < mx , DM can annihilate to the mediator directly and the annihilation cross section is determined 
primarily by the hidden sector coupling gx- 





m\ y \m x J 2 m x m x \ \m x ) 

for the vector and scalar mediators respectively. Meanwhile gj determines how the DM sector couples to the 
SM sector. As for the collider physics, the production of XX occurs through an off-shell mediator; since this 
is a three-body process, the bound is rather weak. Tevatron data requires gf < 0.2 if the mediator couples 



to quarks universally [44 1 



Although gf does not appear to play an important role in the relic density, this coupling controls the 
width (lifetime) of <fi and is relevant for cosmology. The width of the mediator is 



4N f m 4> 2 L f 2m fY fr \ _n» r c m ^ 2 L f 2m f 



f^-YS^V'-W' {T ^ = 2N '^ s T-[^)- (25) 

where the lifetime = T^ 1 . In Section [ill we assumed the DM particles to be in thermal equilibrium with 
the SM thermal bath in the early universe, and in this case the standard freezeout picture and cosmology 
apply. Now, we check the condition for thermalization of the two sectors. If the mediator decay rate is 
larger than the Hubble expansion rate at temperatures T > m^, then the inverse decay processes can keep 
<f> in chemical equilibrium with the SM thermal bath [45]. At these temperatures, the decay rate is given 
by ~ g'jm 2 ^/ (16ttT), where the factor of m^/T accounts for the effect of time dilation. In order for 
the mediator to stay in thermal equilibrium with the SM thermal bath through DM freezeout, we require 
r<£ > H at temperatures T ~ mx- This gives a constraint on gf. 



6 In this paper, we do not consider the intermediate case ~ 2mx , where there is a resonance in the s-channel annihilation 
of XX. 

7 Note that in this case there are also strong bounds on rri0 from neutrino experiments [Hjlh however, we have checked that it 
is still possible to obtain the correct relic density and that the direct detection predictions are unaffected. 
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Figure 3: Lower limit on the mediator mass from combining relic density and DM self-interaction constraints. We 
show the case of a vector mediator; the result for a scalar mediator is similar and is given in Eq. (|32[l . We consider 
DM self-interaction constraints from Bullet cluster observations, elliptical cluster shapes, and elliptical halo shapes. 
The dashed red line indicates the bound on the mass if CMB bounds are also applied, assuming efficiency / ~ 1. 



If g/ is less than the bound given in Eq. (f26|) . the DM sector can have a different temperature from the 
SM sector and the standard freezeout calculation can be modified in a number of ways. We have checked 
that these effects lead to change in the minimum annihilation cross section by less than a factor O(10), 
compared to the results we derived, in Sections iHllIIII Furthermore, the massive mediator is a late-decaying 
particle and in the case where the mediator decays to the SM states, can modify standard nucleosynthesis 
(BBN). There are stringent constraints on the hadronic decay of long-lived particles from the 4 He fraction, 
which requires that the lifetime of the mediator be less than 1CP 2 s 
9q ^ 1-6 x 10~ n ^/l GeV/m^ for a vector mediator, where we take Nj — 3. For leptonic decay modes, we 
take the lifetime of the mediator < 1 s, and obtain a slightly weaker bound, j e >5x 10 _11 -y/l0 MeV/m^, 
for a vector mediator with N'j = 1. 

Finally, we comment on the calculation of the relic density and application of the CMB constraints in the 
light mediator case. When < mx, XX can annihilate to </>(/>, but </> decays to standard model particles 
rapidly compared to the relevant time scales at recombination so that the CMB constraints are unchanged. 
The only difference between a heavy mediator and light mediator with large width is whether there is a 
contribution to the effective degrees of freedom, <?*, from the light mediator. A slightly higher g* in the 
light mediator case gives rise to smaller r^, which in turn weakens the lower bound on (av) from CMB 
constraints. 

In addition, we have neglected the Sommerfeld enhancement effect. As we will discuss in the following 
section, the mediator mass is bounded from below by DM halo shapes; this limits the size of any Sommerfeld 
enhancement. In addition, since (av) rs ira x /m x , for light DM the coupling ax can be much smaller 
and still satisfy the relic density constraint. For the DM masses considered here, we have checked that the 
Sommerfeld enhancement effect is negligible for s-wave and p-wave annihilation processes at both freezeout 
and during recombination, if we take ax and close to their minimum allowed values. 



11 




This leads to a lower bound of 



V. HALO SHAPE CONSTRAINTS ON THE MEDIATOR MASS 



The presence of the light mediator allows for significant DM self-interactions, which can have non-trivial 
effects on DM halo dynamics. A number of astrophysical observations constrain DM self-interactions, for 
example observations of the Bullet Cluster [49} , elliptical galaxy clusters [5(| and elliptical DM halos [HI, [HS . 
Among these, the upper bound on DM self-interaction from the ellipticity of DM halos is the strongest [51] . 
DM self-interactions can erase the velocity anisotropy and lead to spherical DM halos, so the observed 
ellipticity of DM halos constrains the DM self-scattering rate. Because the strength of self-interaction 
increases as the mediator mass decreases, we can use the elliptical halo shape constraint to place a lower 
limit on the mediator mass. Note that in the case of = 0, the ellipticity of the DM halos then places a 
strong upper limit on the hidden sector coupling gx [53| ; it is only possible to obtain the correct relic density 
if m x > 10 3 GeV [Hill^ 8 . 

The effect of DM self-interactions on DM halo shapes can be parametrized by the average rate for DM 
particles to change velocities by 0(1) [H^ : 

T k = J d 3 v 1 d 3 v 2 f(v 1 )f(v 2 ){nxv rcl c7 T ){vl l /v 2 ), (27) 

where nx is the DM density in the DM halo, w rc i = \vi — W2I, and f(v) is the DM velocity distribution in 
the DM halo, for which we take f(v) = e~ v < v ° /(vqs/t:) 2 . ot is the scattering cross section weighted by the 
momentum transfer. 

The form of ot depends on the particle physics nature of DM self-interactions and the relevant momentum 
scales. If the mediator is lighter than the typical momentum transfer in collisions, DM particles interact 
through long-range forces and err depends on velocity. In the opposite limit where the mediator is heavy 
compared to momentum transfer, DM self-interactions are contact interactions and or is independent of v le \. 
In this case, we can take the cxt out of the velocity integrals in Eq. (|27|) and the calculation is straightforward. 
We first will derive the upper bound on the DM self-interaction cross section assuming a contact interaction, 
and then show that this limit applies in deriving the minimum mediator mass. 



We consider the well-studied elliptical galaxy NGC720 [561 |57| . taking our bound from the observed 
ellipticity at a radius of 5 kpc. The DM density profile is fit with local density 4 GeV/cm 3 and radial 
velocity dispersion v r = fg/2 ~ (240 km/s) 2 . We require the average time for DM self-interactions to create 
O(l) changes on DM velocities to be larger than the galaxy lifetime t g ~ 10 10 years i.e. 1 > t g . This 
gives the upper bound 



„ T < 4 , xl0 -^(^)(I^p). (28) 



Other astrophysical constraints have been derived for a /nix, assuming a hard sphere scattering cross 
section a. The bound derived from the Bullet Cluster is (a/mx <2x 10~ 24 cm 2 (mx / GeV)) 49L 5a | ; and 
the bound derived from elliptical galaxy clusters is (a/mx < 10~ 255 cm 2 (mx/GeV)) [50j. 

For the vector and scalar interactions considered here, the force is described by a Yukawa potential 
V(r) = ±axe~ m * r /r. Depending on the mediator, and whether we are in the asymmetric limit, the sign 
may be positive or negative. For the vector case, we have both XX interactions (+) and XX interactions (-) 
unless we are in the asymmetric limit. For the scalar case, the sign is always negative. However, in the limit 
of a contact interaction, the sign of the potential does not matter. We consider i-channel processes in the 
Born approximation. The momentum transfer cross section, defined here by <jt = J d£l*(da/ (Kl*)(l — cos#*), 
is 



27F «2 
°T = -TIP 



ur 



log(l + R 2 " h ' 



l + R 2 



(29) 



8 This limit can be relaxed if the hidden sector is much colder than the visible sector when DM freezes out. In this case, DM 
can achieve the correct relic density with a smaller annihilation cross section [5J|. 
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where /3 — 2axm { j ) /(mxv^ cl ) is the ratio of the potential energy at r ~ m^ 1 to the kinetic energy of the 
particle, and R — rnxv^i/m^ is the ratio of the interaction range to the de Broglie wavelength. For a contact 
interaction R -C 1 , we find that in both cases the momentum- weighted transfer cross section <tt is 

4Tra 2 x m 2 x 

(T T « 4 , (30) 

which is subject to the bound in Eq. (12511 . 

On the other hand, the relic density constraint places a lower bound on the annihilation cross section 
(eru) > 10~ 25 cm 3 /s for light DM and thus on ax- 



for vector and scalar coupling respectively. Note that we assume m^ < mx and take the annihilation cross 
sections in Eq. ([23]). 

Since ax cannot be arbitrarily small, cannot be made arbitrarily small. Combining the bound on ax 
with Eq. (|28p. we obtain a lower bound on the mediator mass: 



m |y > 7 MeV 



(aw) \ 1/4 /m x \3/4 



10- 25 cm 3 /s/ VGeV, 

for the vector and scalar mediator cases, where we take t g = 10 10 years. In Fig. ([3]) we show the lower limit 
on 7710 for the vector case, including the result derived from the more conservative bounds from the Bullet 
cluster or elliptical cluster shapes. We also show the slightly stronger result if we take the CMB constraint 
on the cross section, 9 given in Eq. (fT6|) . 

In deriving the above bound on m^, we have assumed that 3> mxv IC \ and that the Born approx- 
imation is valid. Now we check that the bound given in Eq. (|32p is consistent with these assumptions. 
The condition 3> mxv rc \ is satisfied for 1 MeV < mx < 10 GcV, since from Eq. (|3"2")l we have 
m^/mx ~ 10~ 2 (mx/GeV)~ 1 / 4 but t; rc i ~ 10~ 3 . In this limit the Born approximation is valid if the 
following condition is satisfied 



rV(r)dr 



i) 



mvax 

< 1. (33) 

7770 



From Eq. (pTE]) . we can see v rc \ ax in the DM mass range we are interested in, and thus this condition 
is also satisfied if m^ ^> mxV Te i- We emphasize that we cannot extrapolate the lower mass bound given in 
Eq. ([3"2"1> to mx > 50 GeV because the Born approximation breaks down. For these higher masses, in general 
one has to solve the scattering problem numerically [59[ . In the classical limit where R ^> 1 , there is a fitting 
formula available in [(3(| for the transfer cross section, which has been used to study self-interactions via a 
light mediator for DM masses greater than ~ 100 GeV 45|, 52, 61, 62j |. 



VI. DIRECT DETECTION 

Given the experimental effort needed to detect DM directly, it is important to map out the parameter 
space of direct detection cross sections, subject to the astrophysical and cosmological constraints we have 



9 In the scalar case, annihilation is p-wave suppressed and thus CMB constraints don't apply. 
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discussed. Current experiments are not sensitive to DM-nucleon scattering if the DM mass is below ~1 
GeV because of the energy thresholds. It has been suggested that DM-electron scattering may provide an 
alternative way for the detection of light DM [63[ . We consider DM-nucleon scattering for rax ^ 1 GeV and 
DM-electron scattering for 1 MeV < mx < 1 GeV. 

We compute the range of allowed elastic scattering cross sections within the framework of light DM 
annihilating via hidden sector mediators, assuming mediator couplings to electrons or light quarks. We 
consider both lighter mediators, < mx, and heavier mediators, where we focus on the case 3> mx- 
When rritf, < mx the mediator can be very weakly coupled to the SM, and so the scattering cross sections can 
be much smaller than when 3> mx- However, there is still a lower limit on the cross section coming from 
the lower bounds on the couplings of the mediator to the DM and SM fermions, ax and gj respectively. 
The lower bound on ax is derived from requiring that relic density and CMB constraints are satisfied. 
We consider two possible lower bounds on jj: from requiring the thermalization between the DM and SM 
sectors, or from requiring decay of the mediator before BBN. When 3> mx the lower limit on the cross 
section arises purely from the relic density and CMB constraints. 

Meanwhile, we obtain upper bounds on the electron scattering cross section from the combination of halo 
shape bounds and requiring that the mediator does not significantly affect the electron anomalous magnetic 
moment. Including supernova and beam dump constraints on the dark force coupling 6J| then carves out a 
nontrivial part of the parameter space for electron scattering. 

Fig. ((4]) summarizes our results for the case where the mediator is a vector. We show the possible DM- 
nucleon (left panel) and DM-electron (right panel) scattering cross sections as a function of DM mass. The 
green shaded region is the parameter space for m^ < mx which is allowed by the constraints from the relic 
density, BBN, and DM halo shape constraints; in the electron case we include beam dump and supernova 
cooling constraints. The lighter green area is set by the additional assumption that the mediator has large 
decay width and thus that the two sectors are in thermal equilibrium. In the nucleon scattering case, 
TO0 3> mx is ruled out by CRESST-I and XENON10. In the electron scattering case, the red shaded region 
gives the allowed cross sections for m<f, 3> mx- In the following sections we derive these results and present 
more details. 



A. Nucleon Scattering 

We first consider nucleon scattering in the mass range 1 GeV < mx < 10 GeV, taking universal couplings 
to the light quarks given by g q . The DM-nucleon scattering cross section is given by 

<J n =4a x g 2 n ^, (34) 

where \i n is the WIMP-nucleon reduced mass, and g n = 3g q is the M -nucleon coupling constant. We have 
taken a contact interaction; this is a good approximation because the momentum transfer is much less than 
the minimum mediator mass allowed by the ellipticity of DM halos, as discussed in Section [V] The upper 
bounds here are set by results from direct detection experiments, in particular CRESST-I [(35| and XENON10 

To determine the lower limit on this cross section, we bound ax and g q from below in the case that the 
mediator is lighter than the DM, m<j, < mx- For thermal DM and masses mx > 1 GeV, a lower bound on ax 
is determined primarily by the relic density. As described in Section HTT1 CMB constraints are only important 
in this mass range if </> M decays dominantly to electrons, for which the efficiency factor is / ~ 1. For (j>^ 
coupling primarily to quarks, / ps 0.2 and CMB bounds don't apply above mx ~ 2 GeV. Then the minimum 
annihilation cross section is (av) w -Ka\jm 2 x w 10 _25 cm 3 /s, giving a bound of ax > 5.2 x 10~ 5 (mx/GeV). 
Requiring thermal equilibrium between the hidden and visible sectors, we take the bound on g q in Eq. (|26[) . 
with yJg C R ps 9. Combining the limits above results in a lower bound on the nucleon scattering cross section: 
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Figure 4: (Left) Nucleon scattering through a vector mediator. The green shaded region indicates the allowed 
parameter space of direct detection cross sections. The lighter green region imposes the bound of thermal coupling 
between the two sectors ("large width") while the larger shaded region only requires mediator decay before BBN. 
Also shown is the lower bound for the heavy mediator (m^ ^> mx) case. (Right) Electron scattering through a vector 
mediator, for < mx (green) and m^ 3> mx (red); the intersection of the two regions is shaded brown. We show 
the projected sensitivity of a Ge experiment, taken from [oFj. Beam dump, supernova, and halo shape constraints 
apply here and carve out the region of large <r e at low mx- For more details, see the text. In the lighter green region, 
the condition of thermal equilibrium between the visible and hidden sectors is imposed. 



Since < mx , this quantity is saturated for any mx if we set to its maximum value of m^ ~ mx ■ 
This bound is indicated by the "Large width" line in Fig. Coincidentally, the lower limit here is similar 
to the best achievable sensitivity for WIMP-nucleon scattering if the dominant irreducible background is 



coherent scattering of atmospheric neutrinos off of nuclei 67 69]. However, these studies focused on WIMP 
DM; for light DM, solar neutrinos become much more important and the best achievable sensitivity may be 
several orders of magnitude weaker. 

The lower bound on a n given in Eq. (|35[) is derived by requiring the two sectors be in thermal equilibrium. 
We may relax this assumption, and just demand the mediator decay by nucleosynthesis. This gives g q > 
1.6 x 10 — 11 ^/ 1 GeV/m^, as discussed in Section HV Bl For such g q the two sectors are decoupled through 
freezeout; then the relic density calculation is slightly more complicated and depends on the thermal history 
of the sectors. The change in the relic density then modifies the bound on ax- We have checked that the 
full calculation generally only changes the bound on ax by an 0(1) factor [33], so here we take the bound 
on ax from the large <f> width case for simplicity. In this limit, the lower bound on er„ is given by 



^ K lr ,_ 54 2 / m x \ f GeV\° f fj, 
<r n > 5 x 10 cm x — — - — - 



GSv)uJ (oiiev) (36) 



labeled as "Decay before BBN" in Fig. (g]). 

For reference, we also give the lower bound on the cross section in the case where m^ 3> mx- Here 
DM annihilation occurs directly to SM final states through tp^, with annihilation cross section (av) = 
^ a xgn m x l m %- Since the same combination of parameters enters in both the annihilation cross section and 
the nucleon scattering cross section, we can directly apply the relic density constraint to obtain 
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Figure 5: (Left) Constraints on mediator mass rri0 and coupling to electrons g e for < mx- The shaded region 
is excluded from electron anomalous magnetic moment, beam dump experiments, and supernova cooling 64]. The 
red dashed line shows the g e value used to derive the corresponding red dashed line ("C") in the right plot. (Right) 
Constraints on electron scattering from Fig. [4] The boundaries A, B, and C are discussed in more detail in the text. 



This is the "m^ ^> mx" line in Fig. (j4|). However, this scenario is ruled out by the direct detection limits 
on the cross section. 



B. Electron Scattering 

We consider scattering off electrons for DM in the mass range 1 MeV < mx < 1 GeV. The DM-electron 
scattering cross section is 



ml. 



(38) 



The lower bound on the scattering cross section can be derived in the same way as in the nucleon case, 
taking m$ < mx- Here both CMB and relic density constraints apply, since mx < 1 GeV and the energy 
deposition efficiency f ~ 1 for decay to electrons. We take the bound on the annihilation cross section in 
Eq. (|16[) with c/«l, giving a lower limit on ax- 



ax > 4 x 1CT 7 ( 



mx 



10 MeV 



'hi 



40 GeV 

mx 



(39) 



As in the nucleon case, a lower bound on the DM-electron scattering cross section can be derived by 
assuming that the hidden and visible sectors are in thermal equilibrium. Analogously to Eq. Q35p , we find 



cr e > 3 x 10" 51 cm^ x 



m x 



/10 MeV\ 6 / 

V TTifh I V( 



Me 



2 / /40 GeV\ 
10 MeV J V m 4> J V0.5MeV7 V \ m x )' 



(40) 



where we take y/g c s ~ 3. 

Again, it is possible that the DM sector thermal bath evolves independently from the SM sector and 
in this case we only require the mediator to decay before BBN. From Section IIVB1 we take the bound 
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> 5 x 10 11 ^/10 McV/m^. The minimum scattering cross section is 



* e > 3 x io- cm2 fi^xv (^-)\L (*™\ (4i) 

1 10 MeV/ V TO / V 0.5 MeV/ y V m x / 

If the annihilation goes through a heavier mediator 3> , we derive the strongest lower bound on the 
scattering cross section by applying CMB and relic density constraints: 

^ „ v) 2 A0MeV\ 2 / fi e \ 2 , /40 GeV\ 

a e > 4 x 10" 39 cm 2 ( - - In . (42) 

V m x J V 0.5 MeV/ \ m x ) K ' 

For electron scattering there are no direct experimental bounds on a e . However, for m^ < mx, there 
are bounds on a e arising from indirect constraints, namely halo shape bounds and from searches for new 
light gauge bosons (fjij. The halo shape constraint requires that the self-scattering cross section satisfy 
ax/rax < 4.4 x 10~ 27 cm 2 /GeV with ctt — &iro^m x /m%. If m 4> < mx then constraints on new light gauge 
bosons rule out parts of the (m^, g e ) parameter space; we show beam dump, supernova cooling and electron 
anomalous magnetic moment constraints 10 in Fig. ([5]) (left panel). Here we make use of the convention in 
. where g e = ee, with the kinetic mixing parameter e = eycosOw and e electric charge. The solid line 
(and shaded region) indicates the constraint. 

As a simple application of the constraints discussed above, we derive the upper bound on the cross section 
by rewriting a e : 



4^e / °t ( g, 



^.5xlO-cm 2 (-^-) 2 JH^X. (43) 
~ V0.5McV/ \l m x 

Here we have applied the halo shape constraint and taken (c/e/m^) 2 < 10 _1 e 2 /GeV 2 , arising from measure- 
ments of the electron anomalous magnetic moment (70| . 

To explain more complicated constraints on the (mx,cre) plane from the supernova cooling and beam 
dump experiments for < mx, we show again the allowed parameter space for electron scattering cross 
sections, but highlight boundaries of the constraints by labeling "A", "B", and "C" in the right panel of 
Fig. ([5]). We can map excluded regions on the (m0,g e ) plane to these constraints: 

• Constraint "A": 

For 7710 < mx ;$ 8 MeV, supernova plus beam dump constraints require g e < 1.3 x 10~ 9 . This places 
a stringent upper bound on the cross section, which we derive by taking to its minimum value of 
77i0 = 2m e « 1 MeV, and then setting ax to the maximum value allowed by halo shape constraints: 
a x < 9.5 x 10-VlO McV/m x . This upper bound is then 



a e < 6 X 10" 45 



2 ( y /lOMeV 



Note that the constraint changes somewhat if we also consider m^ < 1 MeV. In this case, su- 
pernova cooling constraints still require g e < 1.3 x 10~ 9 but halo shapes allow for a somewhat 
smaller m^. As a result, the upper bound is slightly weaker if we allow < 1 MeV: a e < 
6 x 10" 44 cm 2 (/Je/0.5 MeV) 2 (10 MeV/mjc) -2 - 



In general there are also constraints from low-energy e+e~ colliders, fixed target experiments, and neutrino experiments [l3j . 
We find these do not significantly affect our results. In the case of kinetic mixing, bounds from measurements of the muon 
anomalous magnetic moment also apply. We do not include them in this paper. 
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• Constraint "B": 

This constraint applies for the large width case. In contrast with constraint A, taking (m^,,g e ) = 
(1 MeV, 1.3 x 10~ 9 ) is in conflict with the condition of thermal equilibrium between the two sectors if 
the DM mass m x > 5 MeV. Furthermore, for m x > 20 MeV, the region (m^ ~ 20 MeV, g e ~ 3 x 10~ 8 ) 
opens up. These competing effects lead to the kink in line B. 

• Constraint "C": 

For mx > 8 MeV, then supernova and beam dump constraints allow a region of larger g e ; for example, 
(m,£ ~ 8 MeV, g e ~ 6 x 10 -4 ) is now allowed. The red dashed lower bound on g e in the left panel of 
Fig. [5] then gives rise to the constraint "C" . The lower bound on the cross section here comes from 
setting ?7i0 ~ mx, applying the red dashed lower bound on g e , and setting ax to its minimum value 
from CMB constraints. 

We make two final notes. First, in the heavy mediator case, the beam dump constraints do not apply and 
the CMB constraints are in general much stronger. As a result, the high cr e , low mx region which is excluded 
in the light mediator case is again allowed indicated by the light red shaded region in Fig. Q. Second, if we 
remove the constraint m^ > 1 MeV, <f> will decay invisibly, and only the supernova constraints are relevant. 
Then a small region of parameter space with g e ~ 1.3 x 10~ 9 and < 1 MeV opens up, as discussed above 
under constraint "A." 

We have verified the bounds discussed above by performing a general scan of the hidden sector parameter 
space. Fig. (jSJ) illustrates our method. We begin by mapping out the parameter space of (ms, g e ) and require 



either large </> width or </> decay before BBN. We combine this with the constraints in [64(, given by the solid 
curve in the top panels of Fig. ([5]). In doing so, we impose the limit 1 MeV < m^ < mx for the case of 
to^ < mx and > 2mx in the case where m^ 3> mx ■ The lower limit of m<p > 1 MeV is imposed in order 
to allow for </> decay to electrons. If the halo shape constraint gives a stronger lower bound on m^,, then 
we take (m^mi^haio < < mx for the < mx case, where (m^) m i n ,haio is minimum mediator mass 
allowed by the halo shape constraint. This generates the sampled points in (m^,, j e ) that we have shown. 
For a fixed (rn^,,g e ), a range of values for ax is allowed, giving rise to a range of allowed scattering cross 
sections. We sample random ax values, subject to the halo shape constraint and the relic density constraint 
as in Eq. (|39p . This then gives a randomly sampled o~ e value, which we indicate by the color of the point in 
Fig. ((5]). For a fixed mx value, because of the range of allowed and ax values, excluded regions in g e 
do not directly map to an excluded region in a e . An excluded region in a e only arises if a sufficiently large 
region of g e is excluded, as shown in the left plot of Fig. ©. We thus verify the possible values of a e in this 
way, imposing all the constraints self-consistently. 



VII. CONCLUSIONS 



Given the unknown nature of DM, it is important to carry out broad-based studies of models of DM. In 
this paper, we have examined constraints on thermal DM with mass 1 MeV < mx < 10 GeV, a mass range 
interesting for numerous phenomenological and theoretical reasons. We considered bounds from cosmology, 
colliders and astrophysics, and derived implications of these constraints on direct detection. 

CMB constraints on DM annihilation present the most serious challenge for light thermal DM, excluding 
symmetric thermal relic DM with s-wave annihilation if mx < 1 — 10 GeV. Two natural ways to evade 
this constraint are to have a DM number asymmetry or velocity suppressed annihilation. In the asymmetric 
case, we found the constraint on the annihilation cross section such that the symmetric component efficiently 
annihilates away; the minimum cross section is larger than the usual thermal relic cross section by a factor 
of a few, depending on the mass. 

Achieving this minimum cross section is difficult if annihilation occurs through a weak scale (or heavier) 
mediator. Collider and direct detection constraints have forced the presence of relatively light mediator 
states in the hidden sector in order to achieve the correct relic abundance and evade the CMB bounds. On 
the other hand, we found that the DM halo shape bounds on DM self-interactions require that the mediator 
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Figure 6: For fixed rax and a mediator with mass < rax, we generate random values of (m^, g e ) allowed by beam 
dump, supernova, a e , and BBN constraints. We show a sample of allowed points in the (m^,, g e ) parameter space; the 
solid curve is extrapolated from the constraints in [r34| . also shown in left panel of Fig. {EJ. For each (m^,<? e ) point, 
we then sample the allowed ax satisfying halo shape and relic density constraints, and compute the corresponding 
elastic scattering cross section a e . The color of the point is determined by a e . (Left) rax = 20 MeV, where the 
minimum mediator mass is ra$ = 1 MeV. (Right) rax = 100 MeV, where the minimum mediator mass m$ > 3 MeV 
is set by halo shape constraints. 



is not too light. We examined constraints from elliptical galaxy NGC720, elliptical clusters, and the Bullet 
Cluster, and derived a lower bound on the mass of the mediator particle. 

We also calculated the range of scattering cross sections allowed within this scenario. Although the lowest 
bound which is cosmologically consistent is well below the reach of any current or envisioned direct detection 
experiments, we showed that several cosmologically interesting benchmarks could be reached. For example, 
in the case of scattering off nucleons, a hidden sector in thermal contact with the SM at T ~ mx can be ruled 
out if an experiment can reach cross sections with a n < 10 _48 cm 2 . In the case of scattering off electrons, 
the scenario where m<f, 3> mx can be probed by direct detection. Beam dump and supernova constraints 
carve out a significant fraction of the available parameter space if < mx- 
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